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Abstract: In this paper, a non- linear ecological mapmay be called as modified Nicholson-Bailey map is considered,which 
becomes chaotic in nature with the increase of the control parameter. As in most of the caseschaos is an unwanted 
phenomenon, so controlling of chaos becomes a necessary part of study. First of all Chan's method is applied on this map 
and the chaotic region is controlled forming periodic trajectories. AgainOGY method is applied on the map to have chaos 
controlled. Lastly, the model has been modified to Chau's form which generates a set of fixed points which have been 
stabilized by OGY method. 
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I. Introduction 

The Nicholson-Bailey model describes the population dynamics of host-parasite (predator-prey) system and is 
described as follows: 

y . — J y p~ a Vn 

yn+i =x„(l-e- ay ") (1.1.1) 

Where x n+1 represents the number of hosts (or prey) at stage n and y n+1 represents number of parasites(or predator) 
at n th stage. It has been observed[21]that the model fulfils the fact that equilibrium state never occurs for predator system in 
nature. Hence a modified version has been introduced by Dutta,T.K.,et al [6] to restrict the unlimited growth of host(or 
prey), which arises in Nicholson Bailey model. The modified model is as follows: 
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y n+1 =x n (l-e- a yn) (1.1.2) 

In fact various modified forms have been discussed in the literature[3,6,ll>14]. 
It has been shown that with the increase of the control parameter, the model (1.1.2) follows a period doubling bifurcation 
route to chaos [6]. As chaos is observed as undesirable part in engineering control. So, a desirable task is to control the 
chaotic region. Since 1990 after the discovery of OGYmethod[22], chaos control problem attracts various researchers. 

We use the periodic pulse method to control chaos which was discussed by N.P. Chau, which restricts the way of 
choosing initial point resulting periodic orbit. Further, OGY method is applied such that the known unstable periodic point is 
made stable with a small disturbance in the control parameter. 

II. Controlling of chaos in the map by Chau's method 

Chau'stheory [2] applied to the map is as follows: 
The model can be written as follows: 

X n+1 = F(X n ), where X n is a vector in R 2 .Let G = KF P , wher K is a diagonal matrix having the diagonal 
elements say k lt /c 2 , and F p is the composition map of F up to p times. If X is a fixed point of G i.e. KF P (X) = X, then the 
fixed point will be stable if the absolute value of the largest eigen value of the Jacobian matrix of G is less than l.The next 
step is to get the value of k lf /c 2 ,such that chaos is controlled. 
For p=l, 



The Jacobian matrix of G is as follows: 
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The fixed points will be stable if \A\ < 1, 
whereA is the eigen value of the Jacobian matrix. 
If X=(x,y) is the fixed point of G, then 

k^Lxe ay x = x , consideringx=£ 0 , we have k x = 



Le -ay-x' 
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(1.2.1) 



Similarly 



k, = 
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x(l-e~ a y) 



(1.2.2) We choose (x,y) such that absolute value of X is less than 1 
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The basin of (x,y) chosen such that it satisfies both (1.2.1) and (1.2.2) is shown graphically as follows: 
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Fig: 1.2.1 Above picture shows the the value of x,y which for particular value of k lt k 2 and p=l becomes a stable fixed point 

for L=3.85 and a=0.1; 
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Fig: 1.2.2: Above figure shows the chaotic region up to 10000 iterations after which chaos control key k lt k 2 is activated to 

make x fixed. 
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Fig: 1.2.3: In the above figure it is clear that y co-ordinate already achieves the fixed value as iteration proceeds. 

Thus as a whole it is clear that once the chaos controlling switch is turned on, the trajectory points (x,y) achieves 
stability to form a fixed point as iteration proceeds. 



If it is desired to obtain a periodic trajectory of period p, we have, 

h = — 



X v-V 
r[ 1 — S: — \-byp-i 
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Where x p .!=first component of F" 1 and y p _i is the second component of f 9 ' . And the Jacobian matrix is given as 

k x 0 \ dx dy 
0 k 2 ) dy^_ dy p 

\dx dy J 
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For p=2; 
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Fig: 1.2.3 Above picture shows the value of x, y which for particular value of k lf k 2 and p=2 becomes a stable fixed point for 

L=3.85anda=0.1: 
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Fig: 1.2.4: Above figure shows the chaotic region up to 10000 iterations after which chaos control key k lt k 2 is activated to 

make x fixed. 
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Fig: 1.2.5: In the above figure it is clear that y co-ordinate already achieves the fixed value as iteration proceeds. 
For p=4, we have 
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Fig: 1.2.6: Above picture shows the the value of x,y which for particular value of k lt k 2 and p=2 becomes a stable fixed 

point for L=3.85 and a=0.1; 

It has been observed that the pool of data as shown in the figure above although serve the purpose of a stable 
periodic point of period two, they are not always able to attract the trajectory towards them as the radius of the basin of 
attraction is small enough. Thus although stable periodic points are created by the kicking method, chaos may not be 
controlled unless the starting point lies in the basin of attraction. 
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Similarly the basin of periodic points of period 8 is created which are stable with the help of a suitable computer program. 

For p=8; 
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Fig: 1.2.7: Above picture shows the the value of x,y which for particular value of k lf k 2 and p=2 becomes a stable fixed point 

for L=3.85 and a=0.1; 

III. Controlling of chaos in the map by OGY method 

The OGYmethod[16,22] applied is as follows: 
Let the two dimensional map (1.1.2) be written as : 

Z n+1 = f(Z n ,L) . Let Z S (L) be an unstable fixed point of equation (1.1. 2) .For values of L near L 0 (say) in a small 
neighborhood ofZ s (L 0 ).The map can be approximated bya linear map given by 

Z n+ i ~ Zs(Lo) = j(Z n - Z S (L 0 )) + C(L - L 0 ) (1.3.1) 

Where J is the Jacobian and C is — , at the point Z S (L 0 ). Assuming that in a small neighbourhood around the fixed 

point 

L — L 0 = —K(Z n — Z s (p 0 )), where K is a constant vector of dimension 2 to be determined. Then the equation (1.3.1) 
becomes 

Z n+ i ~ Zs(Lo) = (J- CK)(Z n - Z S (L 0 )) (1.3.2) 

Using equation (1.3.2) at the parameter value L=3.85, a time series plot is shown below, where after 300 iterations, chaos 
control is switched on. 
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Fig: 1.3. a: Time series graph showing chaos control. 



IV. Applying OGY method on Chau's method 

Let (x*,y*) be a point and let the modified equation of (1.1.2) be 
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Which may be written as h(x,y)=(f(x,y),g(x,y)), where f(x,y)= : - ^ ,g(x,y)= ^^-ay^ clearly h(x*,y*)=(x*,y*).We 
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consider x*=f!(L),y*=f 2 (L), which helps to calculate C = I dg J.Now applying equation (1.3.2) we have 



y 



y n -y 



.(1.4.2) Equation (1.4.2) says that 
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is a fixed point whose stability will be determined the eigen values of J-CK, where K = ^ ^ ^ ^.If X lt X 2 be two 
eigenvalues then k lf k 2 are to be determined in such a way that —1 < 1^ |, |A 2 1 < 1. Calculating J-CK, we have 

J - CK ={ A cl d}^ 
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Let a x = 1 - 2/(L) 2 ,a 2 = /(L)(2/(L)/'(L) + a #'(L)), 
6, = -a/(L),6 2 = /(L)(2/(L)/'(L) + a tf '(L)), 

c l — T777' C 2 — 777^ ' «y V L J — " 



ae -ag{L) g(L ^ 



Then the matrix J-CK= 
Whose eigen values are 

^1,2 
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= 2 ( a i + d i ~ a 2 k i ~ d 2k 2 
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± J(a^ + di - a 2 k! - d 2 k 2 ) 2 + 4 (-%d^ + b^(q - c 2 ki) + (-b^Ci + + b 2 c 2 k 1 )k 2 + sl 2 ^i(^i - d 2 k 2 )) 
Taking x*=y*=L/3,and following the above theory, we have controlled chaos at L=3.85,a=0.1. 
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Fig: 1.4. a representing chaos upto iteration 500 and then control switch is turned on till 800 iteration converging to the stable 

x-coordinate L/3. 
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